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ABSTRACT 

DG Leo is a spectroscopic triple system composed by 3 stars of late-A spectral type, 
one of which was suggested to be a S Scuti star. Seven nights of observations at 
high spectral and high time resolution at the Observatoire de Haute-Provence with 
the ELODIE spectrograph were used to obtain the component spectra by applying a 
Fourier transform spectral disentangling technique. Comparing these with synthetic 
spectra, the stellar fundamental parameters (effective temperature, surface gravity, 
projected rotation velocity and chemical composition) are derived. The inner binary 
consists of two Am components, at least one of which is not yet rotating synchronously 
at the orbital period though the orbit is a circular one. The distant third component 
is confirmed to be a 5 Scuti star with normal chemical composition. 

Key words: Stars: abundances - Stars: fundamental parameters - Binaries: spec- 
troscopic - Stars: variable: 8 Set - Stars: individual: DG Leo 



1 INTRODUCTION 

Different reasons can be invoked to illustrate the fact that 
main-sequence A-type stars occupy a very interesting region 
of the H-R diagram: 

- the transition from radiative to convective energy trans- 
port occurs at effective temperatures between 8500 and 6000 
K (from A5 to F5), a range which encompasses the late A- 
type stars; 

- the classical Cepheid instability strip, when ex- 
tended downward, inters ects the ZAMS (cf. Fig. 8 in 
iRodrfguez fc Bregerlliool at effective temperatures of 8800 
K (blue edge) and 7500 K (red edge) (spectral types from 
A3 to F0); 

- the presence of magnetism is clearly demonstrated from 
the mid-B to early-F spectral types; 

- metal-line abundance anomalies are also frequently de- 
tected among the (non-magnetic) late A-type stars. 

The variety of (atmospheric) phenomena which may affect 
the stars located in that part of the H-R diagram in one or 
the other way is extremely rich: these include several possi- 
ble pulsation driving mechanisms (acting in the <5 Scuti and 
SX Phe, 7 Dor or roAp variable stars) and different processes 
of magnetism, diffusion, rotation and convection which are 
thought to boost or inhibit the presence of chemical pecu- 
liarities in the stellar atmospheres (Ap, Am, p Pup and A 

* This work is based on observations made at the Haute-Provence 
Observatory (France). 



Boo stars). Among the pulsators there exist: the 5 Scuti 
stars which are main-sequence or giant low-amplitude vari- 
able stars pulsating in radial and non-radial acoustic modes 
(p modes) with typical periods of 0.02-0.25 days; the SX 
Phe stars which are high-amp litude S Scuti stars of the old 
disk population jBregerl l2000l) : the 7 Dor stars which are 
cooler and pulsate in non-radial gravity mode s (g modes) 
with typical periods of 1-2 days teregerl lioOO): rapidly os- 
cillating cooler magnetic Ap (roAp) stars that exhibit very 
high-overtone non-radial acoustic pulsations aligned with 
the magnetic axis, inclined relatively to the rotation axis 
('oblique pulsato r' model) an d with typical variability peri- 
ods of 5-16 min dKurt zK OOO'l . Among the chemically pecu- 
liar (CP) stars of interest, we have the magnetic Ap stars 
that are generally variable owing to surface inhomogeneities 
coupled to rotation (model known as the 'oblique rotator'); 
the non-magnetic metallic-lined stars such as the classical 
Am stars which have K-line and metal-line spectral types 
differing by at least 5 spectral subclasses, and the evolved 
Am stars of lumi nosity class I V and III called p Puppis (also 
S Delphini stars. iKurtJ EoOO): and the A Bootis stars which 
show weak met al lines owin g to underabundances of the Fe- 
peak elements (Kurtz 2000). Some of these phenomena may 
or may not be mutually exclusive. For example, enhanced 
metallicity and pulsation can theoretically occur together in 
the cooler and/or evolved part of the instability strip while 
there is mutual exclusion for th e hotter A-stars of the main- 
sequence IITurcotte et al.l2000Tl . The role of multiplicity is an 
additional aspect that has to be considered, as the majority 
of the Am stars are binaries with orbital periods between 1 
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Figure 1. Two hours-averaged cross correlation function ob- 
served from JD2452651.6036 to JD2452651.6954 (i.e.: orbital 
phase ~ 0.98), when all three components are resolved. Each 
correlation peak belongs to one component of DG Leo (Aa, 
Ab or B). 



2 DG LEO 

In this work we present a detailed spectroscopic study of the 
hierarchical triple system DG Leo, which consists of a close 
binary (components Aa and Ab) and one distant compan- 
ion (component B) and whose components were classified 
as late A-ty pe stars. The orbital period of the Aa,b system 
is 4.15 days dPanzieer fc DickenJl967t iFekel fc Bopp||l977tl 
while the orbital period of the visual pair AB w as estimated 
to be roughly 200 years iFekel fc Bopdll977l). Both close 
comp onents have spectral type A8 IV jHoffieit fc Jaschekl 
1982) while the co mposite spectr um was classified as be- 
ing of type F Illn <Cowlevlll976ft and of type A7 III with 
enhanced Sr llCowlev fc Bidelmanlll979Tl . All three compo- 
nents of the system are located in the 5 Scuti instability strip 
and are therefore potential candidates for pulsations. How- 
ever, the claims for short-period oscillations seem mainly to 
concern one component: comp onent B was classified both as 
an u ltra-short-per iod Cepheid ( Eg geiJl979l) and as a S Scuti 
star fclliottll974l) . As a matter of fact, multiple short-period 
oscillations of 5 Scuti type with periodicities of about 2 hrs 
llLampens et aljl2005l) have been detected in the combined 
light curve. 



and 10 days iBudailll996l. Il997l) . Most of these binaries are 
tidally locked systems with synchronized orbital and rota- 
tional periods. Their proj ected equato rial velocities are gen- 
erally below 100 km s" - 1 llKurtd feoOO). Because of this tidal 
braking mechanism, the process of diffusion can act more ef- 
ficiently. In the case of the (magnetic) Ap stars, the braking 
mechanism is magnetic. Their projecte d equatorial veloci- 
ties are usually well below 100 km s" 1 l|Kurtzll200oT) . Thus 
rotation is relevant for explaining the presence of chemical 
peculiarities. An unknown multiplicity status is further ex- 
pected to be the cause of some continuum veiling in several 
A Bootis stars leading to an underestimation of the abun- 
dance of the Fe-peak ele ments as was recently suggested in 
iGerbaldi et"afl l|2003h and lFaraggiana fc Oerbaldil J200?i> . It 
also appears that about 50% of thes e exhibit pulsations of 
S Scuti type llWeiss fc Paunzenll2000l) . Thus complex inter- 
actions between all of these processes exist which are not 
easily disentangled. 

For all these reasons the analysis of the chemical com- 
position of multiple systems having at least one pulsating 
component is an ideal tool to explore in an empirical way the 
interactions that may or may not exist between pulsation, 
diffusion, rotation and multiplicity. In the present paper, we 
are therefore dealing with the chemical analysis of DG Leo 
(HD 85040, HR 3889, HIP 48218, Kui 44), a known mul- 
tiple system with a pulsating component that we describe 
in Sect. |2| Details about the spectroscopic observations we 
performed are given in Sect [3] while the Fourier-transform 
technique that was adopted to disentangle the component 
spectra from the combined ones is described in Sect. 2] Or- 
bital parameters of both spectroscopic and visual orbits are 
discussed in Sect.|^] Fundamental parameters and chemical 
composition of all three components are derived in Sect. HJ 
and Sect. [7| respectively. 



3 OBSERVATIONS AND DATA REDUCTION 

The triple system is spectroscopically unresolved, but the 
Doppler information available in the composite spectrum 
can be used to recover the contributions of the individ- 
ual components. Therefore, extensive spectroscopy was ob- 
tained at the Haute-Pro vence Observatory (O HP) with the 
ELODIE spectrograph fearanne et alJ Il996l) on the 1.93- 
m telescope. Observations were collected during 7 nights in 
2003 (January 3 - 7 & January 11 - 15) and cover about 50 
percent of the close binary's (component A) orbital phase. 
The time exposure was fixed at 360 seconds in order to 
resolve the presumed p-mode pulsations, yet with a good 
signal-to- noise ratio and a high resolution (about 50000). 

The data were automatically reduced order by or der 
with the INTERTACOS pipeline llBaranne et alJll996ft at 
the end of each night. This reduction procedure takes care 
of the order extraction, of the offset and flat-field correc- 
tions and of the wavelength calibration using a thorium ref- 
erence spectrum. The resulting wavelength scale is corrected 
for earth motion afterwards by means of the IRAF soft- 
ware package. In this way, 245 spectra ranging from 3900 
to 6800 A were obtained with a signal-to-noise ratio that 
generally varies between 100 to 180 at 5500 A. 

After each exposure, INTERTACOS provided us also 
with a correlation function (see Fig. computed using a 
template corresponding to a F0 V type star and that ac- 
counts for about 2000 spectral lines. This function allowed us 
to measure accurately the heliocentric radial velocity of DG 
Leo's components, but only at those orbital phases where 
the component structure in the cross-correlation function is 
sufficiently resolved. A sample of the radial velocities mea- 
sured in this way and at these specific phases is listed in 
Table An electronic version of this table is available at 
the Centre de Donnees Stellaires de Strasbourg (CDS). 
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Table 1. Sample of the radial velocities measured on the cross-correlation functions (the 
complete electronic version is available at the CDS). Column 1 gives the Julian date of 
the observations, columns 2, 4 and 6 list the radial velocities of components B, Ab and 
Aa respectively. Errors on radial-velocity measurements (column 3, 5 and 7) are the r.m.s. 
computed by fitting the cross-correlation peak with a Gaussian profile over three different 
intervals. 



JD- 2400000 VB a B VAb cr A b VAa <r Aa 

[km s~ x ] [km s _1 ] [km s _1 ] [km s _1 ] [km s _1 ] [km s _1 ] 



52651.4773 


33.48 


0.52 


-74.64 


0.26 


115.71 


0.56 


52651.4830 


33.16 


0.15 


-76.31 


0.10 


112.27 


0.22 


52651.4886 


35.12 


0.19 


-74.87 


0.05 


112.14 


0.32 


52651.4943 


35.47 


0.18 


-74.60 


0.67 


110.59 


0.37 


52651.5000 


38.14 


0.35 


-74.95 


0.09 


118.18 


0.28 


52651.5057 


38.61 


0.33 


-76.75 


0.06 


117.81 


0.56 


52651.5114 


39.27 


0.22 


-75.58 


0.09 


115.21 


3.10 


52651.5172 


39.53 


0.27 


-75.56 


0.12 


114.48 


0.34 


52651.5229 


39.71 


0.53 


-75.46 


0.10 


114.09 


0.20 


52651.5292 


40.06 


0.68 


-75.66 


0.15 


110.95 


0.74 



Table 2. Sample of radial velocities derived with KOREL (the complete electronic version is 
available at the CDS). Column 1 gives the Julian date of the observations, columns 2, 4 and 
6 list the radial velocities of components B, Ab and Aa respectively. The deviations for each 
measurement relatively to the orbital solution, or (O— C)s, are given in columns 3, 5 and 7. 



JD-2400000 VB (O-C) VAb (O-C) VAa (O-C) 

[km s -1 } [km s _1 ] [km s _1 ] [km s _1 ] [km s -1 ] [km s _1 ] 



52651.4773 


37.60 


-0.87 


-75.01 


-0.52 


112.56 


-2.30 


52651.4830 


36.71 


-1.76 


-76.01 


-1.23 


113.49 


-1.66 


52651.4886 


36.90 


-1.57 


-75.90 


-0.84 


114.02 


-1.41 


52651.4943 


37.34 


-1.14 


-76.75 


-1.40 


115.45 


-0.25 


52651.5000 


37.62 


-0.85 


-76.13 


-0.51 


114.74 


-1.24 


52651.5057 


38.28 


-0.19 


-76.72 


-0.83 


115.49 


-0.75 


52651.5114 


38.47 


0.00 


-77.49 


-1.35 


113.88 


-2.62 


52651.5172 


38.83 


0.36 


-76.73 


-0.34 


115.73 


-1.02 


52651.5229 


39.14 


0.67 


-77.07 


-0.43 


115.53 


-1.47 


52651.5292 


38.47 


0.00 


-77.20 


-0.30 


114.75 


-2.51 



4 SPECTRAL DISENTANGLING 

The composite spectra of DG Leo are complex, and a care- 
ful quantitative analysis is needed. To extract the individ- 
ual contributions of the three components, we adopted the 
spectral disentangling technique. This technique determines 
the contributions of the components to the composite spec- 
tra and the orbital parameters in a self-consistent way. The 
feasibility of the disentangl i ng me thod in practice was first 
proved bv ISimon fc Sturml lll994j) . They reconstructed the 
spectrum of the primary component of V453 Cygni from 
composite spectra at different orbital phases out of eclipse 
and showed that it matched perfectly the one observed dur- 
ing total eclipse in the secondary minimum. Since then, dif- 
ferent implementations of the method have been success- 
fully applied to binaries wit h components ranging fr om O- 
type (ISturm fc Simonlll99l to F-type components jGriffinl 
2002). Advantages and disadv antages of se veral implemen- 
tations have been discussed bv llliiid i2004t) . 

We used the technique as developed by lHadraval 



il995l and references therein) and applied in the korel com- 
puter code (Release 21.3.99). Although the method no longer 
requires the intermediate step of deriving radial velocities, 
KOREL cross- correlates the input spectra with the disentan- 
gled component spectra afterwards to provide relative radial 
velocities for comparison with historical data. Note that the 
systemic velocity must be determined independently, as the 
concept of spectral lines does not enter in the purely math- 
ematical disentangling procedure. A basic assumption on 
which the disentangling is carried out is that the shape of the 
line profiles of the binary components should not vary with 
time. However, the technique has been found to be applica- 
ble in the presence of complex variations of l i ne shapes un- 
der c ertain conditions llHarmanec et alll2004t ll)e Cat et alJ 
120041) . Helpful conditions are: a very different time-scale for 
line variability and orbital revolution; many spectra and line 
variability not phase-locked with the orbital motion; ampli- 
tude of line variability small with respect to the weakest 
contributing component. Under these conditions the varia- 
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Figure 2. Radial velocities obtained with KOREL after the dis- 
entangling procedure, and reported in the heliocentric frame- 
work. 

tions simply act as additional high-frequency noise. We are 
in a relatively favourable case in all aspects and especially 
because the three components each contribute roughly one 
third to the composite spectrum while line variability will 
turn out to be moderate and limited to only one component. 
All 245 spectra obtained at OHP were used in the disentan- 
gling procedure. 

4.1 Component spectra and luminosity ratios 

When the components' relative fluxes are known or when 
they v ary substantially with t ime (e.g. when eclipses are ob- 
served, |HeMberg^^i]|20o3) , the method does allow to ac- 
count straightforwardly for the veiling owing to multiplicity 
which artificially weakens all the spectral lines. Illiiic et all 
( 2004) show that in the case of time-independent component 
fluxes, the spectral disentangling can be performed assuming 
equal fluxes and the resulting spectra can be re-normalised 
afterwards. Therefore, we first assumed that each compo- 
nent contributes equally to the total luminosity. Then, to 
correct the resulting component spectra, we estimated the 
relative luminosity at A3930 A with the assumption that the 
Ca II K line depth is saturated fSect. lS"3*l . We finally applied 
this ratio to the whole investigated wavelength range in view 
of the similar flux distribution of the 3 stars (Table |7J. As 
a result, the luminosity ratios and the estimated individual 
V magnitudes we derived are given in Table |3] The implied 
difference between the magnitudes of components A and B 
(AV = 0.59 ±0.09) is consistent with the differential m agni- 
tude measured by HIPPARCOS (AHp = 0.69 ± 0.08; ESA 
119971) . Examples of the resulting components' spectra can 
be found in Fig. |S| The complete disentangled spectra are 
available at the CDS as ps-format files. 

4.2 Radial Velocities 

As noted earlier, KOREL cross-correlates the input composite 
spectra with the output disentangled component spectra af- 
terwards in order to provide radial velocities with respect to 
the centre of mass of the system. To estimate the systemic 
velocity of the DG Leo system, the radial velocities deduced 
with KOREL were compared to those measured on a subset 



Table 3. Luminosity ratios (rx,) and compo- 
nent magnitudes (Vi) derived from the depth 
of the Can K line. Note that errors are not 
uncorrelated, since the combined magnitude is 
known accurately and the sum of the fractional 
luminosities is exactly unity. 



Component 
















(%) 




(mag) 


DG Leo Aa 


32 


± 2 


7.31 


± 


0.15 


DG Leo Ab 


31 


± 2 


7.35 


± 


0.15 


DG Leo B 


37 


± 2 


7.16 


± 


0.14 



of spectra with INTERTACOS (7 = 26.54 ± 0.5 km s _1 ). 
A sample of these KOREL-based heliocentric radial velocities 
is given in Table [5] The complete data set is electronically 
available at the CDS and plotted in Fig.[2]against the orbital 
phase of the close binary (Aab system). 

A detailed comparison of both types of radial-velocity 
measurements (Tables0and|5J in the common subset shows 
no significant offset (at the 1 km s~ level) for any compo- 
nent, but a different scatter for the three stars: the largest 
scatter (differences up to 10 km s^ 1 ) is observed for the pul- 
sating (Sect. 8.6) B component, while the smallest scatter 
(differences up to 3 km s _1 ) applies to the Ab component. 
The wider lines of Aa produce differences up to 6 km s -1 . 
The more realistic correlation masks from the KOREL proce- 
dure undoubtedly produce better radial-velocity estimates 
than the Gaussian component fitting of the INTERTACOS 
CCFs (in the sense of lower residuals relative to the orbital 
solution), especially for the partly overlapping Aa and B 
components. Therefore, the radial velocities used in the fur- 
ther analysis are those derived for all 245 spectra with the 
component spectra as correlation masks. 



5 ANALYSIS OF THE ORBITS 

5.1 Spectroscopic orbit of the close binary Aab 

The orbital parameters were derived applying KOREL to all 
245 spectra in 5 spectral regions. The regions were selected 
for showing the largest flux gradients. These wavelength do- 
mains were basically those having simultaneously the great- 
est line density and the fewest blends. During the fitting 
procedure, 5 parameters were finally used to describe the 
orbit of the close binary system (A) and were defined as free 
parameters, while the longitude of periastron was kept fixed 
(i.e.: uj = 0°). The final solution is: 

P = 4.146751 ± 0.000005 days 

T = 2452639.259 ± 0.001 JD 

e = 0.0000 ± 0.0004 

K Aa = 100.72 ± 0.18 km s _1 

q = 1.000 ± 0.001 

where P is the orbital period, To is the time of nodal passage, 
e is the eccentricity, Kah is the amplitude in radial velocity of 
component Aa and q is the mass ratio. Notice that, as the or- 
bit is circular, the time of nodal passage was adopted. Since 
KOREL does not give error estimates, uncertainties were esti- 
mated from the covariance matrix using the radial- velocities 
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as input. This is a conservative way of estimating uncer- 
tainties, since it involves the derived radial velocities, that 
are vulnerable to cross-correlation bias in a composite spec- 
trum. Actually, plotting the radial-velocity residuals (listed 
in Table |5J against phase indicates correlated errors at the 
1 km s _1 level, with strings of characteristically 6 consecu- 
tive (O— C) having the same sign. This is taken into account 
explicitly in the uncert ainty on the velocity amplitu de by 
the method proposed bv lSchwarzenberg-Czernvl lll998h . Fur- 
thermore, the accuracy of the orbital period was obtained 
by taking into account the uniquely determi ned number of 
cycles between our observations and those of iFekel fc BqppI 
( 1977). It is further worth noting that the period of the vi- 
sual binary (AB) is so long compared to the time coverage of 
our observations that its orbit and motion can be neglected. 



5.2 Astrometric— spectroscopic orbit of the wide 
binary AB 

Altho ugh micromet ric observations of DG Leo exist since 
1935 jKuiperlll935l) . no orbit is known yet. Fek el fc BoppI 
Jl977h roughly estimated an orbital period of 180 ± 50 yrs 
based on calibrations of spectral type versus mass and ab- 
solute magnitude. Notice that in their study a parallax of 
0.007"was adopted which is in full concordance with the 
Hipparcos parallax (6.34 ± 0.94 mas). In order to derive 
the orbital parameters, we requested all the observations 
of DG Leo AB f rom the Washingto n Double Star Obser- 
vations Catalog jMason et al-lEoOlh and obtained 44 mi- 
crometric and 38 speckle-interferometric measurements (see 
Table |SJ. We further received one unpublished speckle ob- 
servation (Hartkopf, W., 2004, private comm.). Since 1997 
the systematic observations were discontinued owing to the 
difficulty of access to sufficiently large telescopes for this 
type of object. Unfortunately, the astrometric data in them- 
selves are insufficient to compute a reliable visual orbit, as 
practically no change in position angle was detected dur- 
ing a time interval of more than 40 years. To complete the 
data we therefore included the components' radial velocities 
(see Table^J for a combined astrometric-spectroscopic orbit 
analysis. Use wa s made of the software package developed by 
iPourbaixl jl998l). It is based on the pr inciples of simulated 
annealing ("SA. ^Metropolis et alj ll*953) for the global explo- 
ration and minimization in the parameter space followed by 
a local least-squares minimization (following the scheme of 
Broyden-Fletcher-Goldstrab-Shanno). Table |S| lists the sta- 
ble and well-defined orbital parameters corresponding to one 
of the best solutions in the sense of least-squared residuals 
after 1500 iterations with SA and exploring the period range 
70-200 years. Fig.[3]illustrates this orbit, both for the visual 
(upper panel) and the spectroscopic (lower panel) datasets 
(cf. Table The accuracy and the confidence level of this 
solution are further discussed in Sect. 



6 ATMOSPHERIC FUNDAMENTAL 
PARAMETERS 

The disentangled component spectra provided by korel al- 
low the use of well known techniques of spectral analysis 
that lead to the determination of physical quantities (i.e.: 



Table 4. Radia l vel o cities of component s A and B. 
FB77: IFekel fc BopeI Il977ll : RS91: iRosvick fc Scarfe j 
EB2002: ELODIE observations made bv lOblalJ 
ELODIE observations performed by 
20031): PM2004: ELODIE observations ob- 



1991 

20021) : PN20 03 
iNorthl 



|200 



taincd by Mathias (2004). 



JD 


Va 


Vb 


Ref. 




(km s" 1 ) 


(km s — 1 ) 




2440995.9042 


26.9 ± 3 


28.1 ± 5 


FB77 


2442113.9789 


25.1 ± 2 


27.9 ± 4 


FB77 


2442470.8851 


27.2 ± 3 


29.5 ± 3 


FB77 


2442856.8086 


29.3 ± 2 


24.3 ± 4 


FB77 


2447815.7210 


29.9 ± 2 


22.2 ± 2 


RS91 


2452300.6250 


19.3 ± 2 


41.6 ± 4 


EB2002 


2452649.0512 


20.1 ± 1 


38.5 ± 4 


this work 


2452771.3452 


21.1 ± 4 


43.0 ± 5 


PN2003 


2453026.6498 


22.0 ± 4 


32.0 ± 5 


PM2004 



Table 5. Sample of micrometric (9, p) and speckle 
(x, y). The data are followed by the residuals (O— C). 
A full version of this table is available at the CDS. 



MICROMETRIC DATA 



Epoch 


en 


(0-C)6» 


Pi") 


(O-C)p 


1935.4400 


218.360 


+3.548 


0.335 


-0.005 


1936.8700 


218.118 


-5.715 


0.339 


-0.019 


1937.1100 


218.078 


+11.524 


0.339 


-0.039 


1937.2800 


218.049 


+8.852 


0.340 


-0.040 


1938.1600 


217.904 


+0.195 


0.342 


-0.042 


SPECKLE DATA 


Epoch 


x(") 


(O-C)x 


y(") 


(O-C)y 


1940.1200 


-0.274 


-0.010 - 


-0.211 


+0.008 


1976.8606 


-0.249 


-0.001 - 


-0.152 


-0.004 


1977.1802 


-0.247 


+0.011 - 


-0.151 


-0.032 


1977.3331 


-0.247 


+0.002 - 


-0.150 


-0.006 


1978.3161 


-0.242 


-0.015 - 


-0.146 


+0.001 



Table 6. Values and standard deviations of the 
constrained parameters of the combined astrometric- 
spectroscopic orbital solution of DG Leo AB. 

Orbital parameter Value Std. dev. 



«(") 


0.191 


0.015 


i P) 


117. 


13. 


w n 


341. 


13. 


np) 


27.2 


3.0 


e 


0.946 


0.043 


p (yr) 


102.3 


9.5 


T (Besselian year) 


1999.59 


0.41 


V (km s" 1 ) 


+26.97 


0.49 


K ~ M A +M B 


0.374 


0.046 
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Table 7. Derived fundamental parameters 



-0.2 -0.1 
A Y (arcsec) 





Orbital solution 1 (e - 0.95) 

Orbital solution 2 (e - 0.99) 




1970 1980 1990 2000 

Bessel year 

Figure 3. Upper panel: Observations and computed relative 
visual orbits: the best solution is represented by the thick solid 
curve while the other represents an orbit with an eccentricity very 
close to 1. The micrometric and speckle data obtained from 1935 
to 1997 are represented by open and filled circles respectively. 
Lower panel: Observed (comp A: open triangles — comp B: 
filled triangles) and predicted (comp A: solid curve and comp B: 
dashes) radial-velocity variations for the two solutions (thick and 
thin curves) shown in the upper panel (see discussion in Sect, 18.21 . 



T c ff, log g and Tasini) accurately describing the stellar at- 
mospheres. For DG Leo's components, these atmospheric 
fundamental parameters were derived in three consecutive 
steps which are chronologically described in the following 
subsections and the results of which are listed in Table 



6.1 Projected rotation velocity 

A Fourier analysis of several unblended line profiles was 
used to determine the projected rotation velocity (V sin i). 
The adopted proce dure that us es a FFT FORTRAN sub- 
routine pr ovided b y iGraj «1992f) is similar to the one de- 
scr ibed bv | Carro l|j[l933l) and largely applied for example 
by iRover et alJ i2002l) . It is based on the location of the 
first minimum of the Fourier transform that varies as the 
inverse of V sin i. Twelve isolated spectral lines were there- 
fore chosen in the disentangled spectra (see Sect. 4). The 
projected rotation velocities thus derived are given in Ta- 
ble Q The listed error is equal to the root mean square of 
the measurements made on each of the 12 unblended lines. 
Our Vsm i measurements confirm the general trend that is 



Star 



Aa 



Ab 



T cff (K) 
log 9 

V sin i (km s _1 ) 
€turb (km s" 1 ) 
M H r (M q ) 
Rhr (R©) 



7470±220 
3.8±0.14 
42±2 
2.3±0.5 
2.0±0.2 

2.96±0.20 



7390±220 
3.8±0.14 
28±2 
2.3±0.5 
2.0±0.2 

2.94±0.20 



7590±220 
3.8±0.12 
31±3 
2.5±0.5 
2.1±0.2 

2.99±0.20 



observed when the widths of the cross-correlation functions 
(Fig. for each component are compared to each other: 
DG Leo Aa is apparently rotating faster than DG Leo Ab. 
This unexpected result is discussed in Sect. 8.5. 



6.2 Effective temperature 

In the atmospheres of mid- to late-A-type stars, hydro- 
gen exists mostly in a neutral form. Hydrogen lines are 
therefore rather insensitive to surface gravity and can be 
used to derive the effective tem perature (see for exam- 
ple: ISmallev fc Dworetskvl Il99^). Consequently, we esti- 
mated the effective temperature of DG Leo's components 
by fitting the individual Ha and H7 lines with theoretical 
line profiles (Fig.^J. Fully line-blankete d model atmosp heres 
were computed in LTE with ATLAS 9 (lKurucdll993al) cor- 
rected to account for the comments made bv lCastelh et alJ 
1119971) . From these models, a flu x grid was obtained by use 
of the synspec computer code llHubenv fc Landll995l see 
references therein) enabling the IRSCT and IOPHMI opac- 
ity flags, respectively, to account for the Rayleigh scattering 
and for the H~ ions. A least-squares method and the Ml- 
nuit minimization package of CERN were finally used to 
fit the observed line profiles. During the procedure, the sur- 
face gravity was kept fixed to 4. (cgs) for all three compo- 
nents. Error bars were computed from the largest deviation 
obtained by performing the fit several times and assuming 
different start values for the effective temperature. 



6.3 Surface gravity 

Because in late-A-type stars the hydrogen lines are in- 
sensitive to a change in surfa ce gravity, we fol lowed 
the same procedure as used by iNorth et alJ <ll997l) and 
lErspamer fc North! (120031) for stars within a distance of 150 
pc to derive the surface gravity of the components. The lu- 
minosity was calculated from the HIPPARCOS parallax of the 
triple star and the components' V magnitude s (see T able Et. 
The bolometric correction was adopted from iFloweJ |l996). 
Making use of the bolometric magnitude and the previously 
derived effective temperature, we then obtained the mass 
and the radius of each component th rough interpola t ion in 
the theoretical evolutionary tracks of ISchaller et alJ dl992l) 
(for Z=0.001 and 0.020). In this way the surface gravities 
of the components were found to be much alike: the values 
are listed in Table [7] together with the interpolated masses 
(Mhr) and radii (Rhr). Their accuracy was estimated from 
the errors affecting the HIPPARCOS parallax (^=15%) and 
the effective temperature. 
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7 ABUNDANCE ANALYSIS 

The disentangled spectra were used to estimate the compo- 
nents' atmospheric chemical composition. We therefore sys- 
tematically applied korel, order by order and from 4900 
A to 6100 A, keeping the orbital parameters and luminos- 
ity ratios fixed to the values in Tableland Sect. IS"T1 The 
abundance determination was then carried out by fitting 
theoretical line profiles to the observed components' spectra 
(see Fig. 

7.1 Atomic data 

Oscillator strengths, energy levels and damping parameters 
(including Stark, van der Waals and natural broadening) 
used during the fitting procedu re were basically th ose com- 
pil ed in the VALD-2 databas e jKupka et al] ll999) updated 
bv lErspamer fc North! i2002T) . For elements having atomic 
numbers greater than 30, the line list was compl eted by 
atomic data extracted from Kurucz CD-ROM N23 iKurucd 
Il993bl) and from the NIST A tomic Spectra Databas e (see 
the gfVIS.dat file created bv iHubenv fc Land Eoblh . Sen 
an d Mnu were treate d with oscillator strengths proposed 
bv lNissen etafl (|2000D . 

7.2 Chemical composition 

The microturbulent velocity was supposed to be constant 
with optical depth, ft was derived simultaneously with the 
iron abundance by fitting the observations, from 4950 A to 
6100 A, with theoretical spectra. Except for iron, during this 
procedure the chemical composition was kept fixed at the 
solar values. The resulting microturbulent velocities, £turt>, 
are given m Table □ and were further used to derive the 
abundances of the other atomic species. 

A description of the model atmospheres and of the 
computer code we used can be found in Sect. 16.21 For 
the same line lists and model atmospheres, there is gen- 
erally a good agreement between the synthetic spectra pro- 
vided by SYNSPEC46 and those computed by other programs 
such as spectrum wh ich was written by R.O. Gray (e.g.: 
iGrav fc Corbalhll994D . However, Sc II line profiles were sys- 
tematically much stronger in SYNSPEC's spectra, providing 
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Figure 5. The abundance patterns we observed for the compo- 
nents of DG Leo are com pared to those determine d in the Sun 
jGrevesse fc Sauvall99sl) and in the Am star 32 Aqr llKocer et alj 
Il993t) . 



scandium abundances one order of magnitude smaller than 
expected for our reference star. The use of Kurucz's PFSAHA 
subroutine, to compute the partition functions of scandium, 
definitely resolved this incoherence. It is worth noting that 
this problem does not occur with the older versions of SYN- 
SPEC. 

The final results of the chemical analysis are summa- 
rized in Table |H|and in Fig. [^J An example of the agreement 
obtained between the fitted synthetic spectra and the dis- 
entangled ones is shown in Fig. HJ Logarithmic abundance 
values are given relatively to hydrogen for the Sun (Col. 6) 
as well as for each component of DG Leo (Col. 2, 3 and 
4). So lar abundances were taken from iGrevesse fc Sauvall 
(1998). The number of transitions considered to derive the 
atmospheric chemical composition of DG Leo's components 
is noted in Col. 5. Abundance error bars, Ae, were computed 
in the framework of non-correlated error sources with the 
following relation: 



8 Y. Fremat, P. Lampens & H. Hensberge 



Table 8. Components' chemical compositions relative to hydro- 
gen 



spectroscopic orbital parameters listed in Sect 5.1, one may 
further derive: 





e(Aa) 


DG Leo 

e(Ab) 


6(B) 




Sun 


Na 


-5.47 ± 0.13 


-5.52 ± 0.13 


-5.63 ± 0.13 


1 


-5.67 


Mg 


-4.42 ± 0.17 


-4.64 ± 0.17 


-4.27 ± 0.17 


2 


-4.42 


Si 


-4.48 ± 0.11 


-4.51 ± 0.06 


-4.48 ± 0.05 


3 


-4.45 


Ca 


-5.65 ± 0.15 


-5.84 ± 0.17 


-5.57 ± 0.13 


4 


-5.64 


Sc 


-9.29 ± 0.16 


-9.35 ± 0.16 


-8.87 ± 0.16 


2 


-8.83 


Ti 


-6.74 ± 0.09 


-6.87 ± 0.09 


-6.80 ± 0.08 


5 


-6.82 


Cr 


-6.03 ± 0.12 


-6.12 ± 0.12 


-6.31 ± 0.12 


8 


-6.33 


Mn 


-6.66 ± 0.13 


-6.60 ± 0.13 




2 


-6.61 


Fe 


-4.32 ± 0.07 


-4.33 ± 0.08 


-4.54 ± 0.08 


30 


-4.50 


Ni 


-5.40 ± 0.09 


-5.50 ± 0.11 


-5.63 ± 0.11 


7 


-5.75 


Y 


-8.78 ± 0.14 


-8.67 ± 0.14 


-9.74 ± 0.16 


5 


-9.78 


Zr 


-8.49 ± 0.16 


-8.39 ± 0.16 




1 


-9.39 


Nd 


-9.50 ± 0.26 


-9.48 ± 0.26 




4 


-10.50 



Table 9. Abundance deviations (Ae) caused by the upper and 
lower limits of T e fj , log g and microturbulent velocity (£turb ) • 



Elem. 


Ae(Teff) 


Ae(log i 


l) Ae(£ turb ) 


Na 


0.03 


0.02 


0.05 


Mg 


0.10 


0.04 


0.07 


Si 


0.02 


0.02 


0.00 


Ca 


0.03 


0.01 


0.06 


Sc 


0.01 


0.06 


0.13 


Ti 


0.02 


0.05 


0.05 


Cr 


0.06 


0.01 


0.03 


Mn 


0.05 


0.00 


0.04 


Fe 


0.01 


0.00 


0.06 


Ni 


0.02 


0.03 


0.05 


Y 


0.02 


0.05 


0.09 


Zr 


0.02 


0.11 


0.05 


Nd 


0.08 


0.07 


0.05 



(a Aa + a A b) siniA = 
M Aa sin 3 i A = 



16.5 ± 0.02 R 
1.760 ±0.007 M Q 



= M 



Ab sm i A 



(2) 
(3) 
(4) 



where i A stands for the inclination of the close binary and 
where we introduced the values of K, q, P and e as well as 
their a ssociated error bars. Since no eclipses were observed 
iLampens et alj|2005l) . the derived orbital parameters con- 
strain the orbital inclination angle which will depend on the 
stellar radius and thus on the surface gravity. From the con- 
dition R Aa + R A b < (a Aa + a A b) cosi (no eclipses), we have 
that: 



tan i A < 



("Ai+'Ab) 



RA a +RAb 

For R Aa = R Ab = R and if we use R" 1 cx (g/M) 1/2 in 
relations J2J and ©, we obtain: 



. . 2 • 
sin i a cos i A 



> 
> 



Msin 3 i A /M Q g s 
((a Aa + a Ab )/2) 2 g 

0.02586 — 



Numerically for log 5 = 3.8 ± 0.14 (see Table |7J, we find: 



(•>) 

(6) 

.: 
(7) 



implying an upper limit for the orbital inclination of 73° 
from the absence of eclipses. Comparing the stellar masses 
derived from the evolution tracks (see Mhr in Table|7J, it is 
clear that the orbital inclination cannot be much lower, as 
the lower mass limit deduced from (J7J and @ gives: 



M Aa = M Ab > 



2.11°;? M Q . 



(8) 



a 2 = ^cai. + Ae(Teff) 2 + Ae (logg) 2 + Ae(£ turb ) 2 (l) 

n — 1 

where n is the number of lines that we used and where <r sca t , 
Ae(Teff), Ae(log g) and Ae(£ t m-b) are the errors related re- 
spectively to the r.m.s. scatter of the derived abundance, to 
the effective temperature, surface gravity and microturbu- 
lent velocity accuracy. The uncertainties Ae on the abun- 
dance introduced by the fundamental parameters are given 
in Table [5] They were computed relatively to a reference 
atmosphere model defined by the following fundamental pa- 
rameters: T cf f = 7500 K, log g — 4.0, £ tU rb = 2.5 km s" 1 and 
solar chemical composition. When only one line is available 
to determine the abundance, the computation of o" scat is 
based on the oscillator strength's relative error (fixed at 25% 
when not found in the NIST database). 

8 DISCUSSION 

8.1 Spectroscopic orbit of the close binary Aab 

The parameters we derived for the spectroscopic binary 
are in good agreement with the previously pub lished val- 
ues jFekel fe Bopdll977t iRosvick fe Scarfdll99ll) . From the 



8.2 Astrometric— spectroscopic orbit of the wide 
binary AB 

Recent and old astrometric data were combined with radial 
velocities obtained at an epoch close to periastron passage. 
Using different starting points corresponding to the 100 low- 
est values of the objective function found with SA, we ended 
up with the solution proposed in Sect. 5.2 in 80% of the cases 
after a fine tuning of the minimization process. In only 3% of 
the cases, SA converged successfully to solutions with some 
parameter out of the 1.5-cr uncertainties given in Table |S| 
Although the residuals are realistic given the measurement 
errors, the orbit is not yet definitive. Some parameters re- 
main unconstrained because we miss radial velocities near 
periastron in a very eccentric orbit. At present, we conclude 
that: 

- the orbit al period is notably shorter than previously 
estimated by iFekel fe Bopd lll977l) , with a most probable 
value close to 100 years (and an uncertainty of 10%); 

- the orbit is definitely very eccentric, but unfortunately 
we do not have radial velocities near periastron. Any spectra 
taken in the interval 1998-2002 would be of immense value 
if existing; 
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Figure 6. Comparison between the vertically shifted normalized component spectra (thick line) and the fitted synthetic spectra resulting 
from the chemical abundance analysis (thin line). Identifications indicate the transitions contributing the most to the spectral lines. 
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- the fractional mass is constrained reasonably well with a 
most probable value k = 0.37, and consistent with the value 
derived from the fundamental atmospheric parameters (k — 
0.34 ± 0.07). The constraint on the orbital k already induces 
a useful lower limit to the mass of component B relative to 
each of the equal-mass components of A: 



Me 
M a 



> 1.19" 



(9) 



In other words, component B is very probably at least as 
massive as each of the close binary components; 

- the other constrained parameters (Table |SJ include the 
orbital inclination (range 110-130 degrees) and those related 
to the astrometric data; 

- ambiguity remains with several parameters related to 
the unconstrained radial- velocity amplitude (see the two dif- 
ferent solutions shown in Fig. 3). As a consequence, we do 
not know the true semi-axis major (in km s _1 ), which would 
provide a dynamical parallax and absolute masses. 

Hence, we can use the mass of component A obtained 
from the close spectroscopic orbit. From Table |7| and eq. 8, 
we have Ma = 4.2 ± 0.2 Mq. Then the data of Table EJ 
together with standard-error propagation laws, lead to: 



Mb 
Ka 
Kb 
a A + a B 

TV 



= 2.5±0.6M Q 

= 12.3 ± 5.8 km s" 

= 20.6 ± 8.7 km s" 

= 41.2 ± 5.8 AU 

= 4.6 ± 1.0 mas, 



where the dynamical parallax is consistent with the Hippar- 
cos parallax at the 1.3— tr level. 

A next step would be to repeat the SA procedure limit- 
ing the exploration of the parameter space to realistic values 
of Ma- However, significant progress in the direct determi- 
nation of the individual masses with a sufficiently high accu- 
racy could be achieved, either by including possibly existing 
radial velocities near periastron passage - we therefore invite 
anyone who might have spectra or radial velocities of DG 
Leo during the interval 1998-2002 (see Fig. |SJ to contact us 
- or by gathering new speckle-interferometric observations 
in the following years. 



8.3 Ca abundance and luminosity ratios 

As noted in Sect. 01 the luminosity ratios of the three compo- 
nents were derived under the assumption that their Ca II K 
line-depths are saturated and identical. To test the validity 
range of that assumption, we compare in Fig. 0the theoret- 
ical spectra of all three components of DG Leo from 3920 to 
3950 A. Owing to the close values of the fundamental pa- 
rameters and to the fact that the calcium underabundance 
observed in the close binary components is smaller than 0.20 
dex, the line depths are very similar. As a matter of fact, the 
slightly smaller depth in the spectrum of the faster rotating 
component Aa might have introduced a flux zero-point error 
smaller than 1.5 percent, which does not significantly affect 
the abundance analysis. 








□G Leo Aa 
DG Lea Ab 
DG Leo B 
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3930 3940 
Wavelength (A) 



3950 



Figure 7. Computed Ca II K line-profile for all three components 
of DG Leo. 
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Figure 8. The spectrum of the chemically normal star HD 33959 
(top) and the ratio of the Am and DG Leo spectra (components 
Ab and B only) to the top one are plotted. 
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Figure 9. The spectra of DG Leo Aa (line) and DG Leo Ab 
(dots) are compared. The spectrum of the latter component was 
convolved with a Gaussian function in order to increase artificially 
its projected rotation velocity. 



8.4 Comparison with reference stars 

If we discount the larger effects of Doppler broadening for 
DG Leo Aa, both components of the close binary have very 
similar spectra. Their nature is however different from that 
of component B and their iron-peak lines are generally sig- 
nificantly stronger. To visualize the general trend that can 
be noted in the chemical abundances we derived, their re- 
spective spectra have been compared to the spectra of the 
A9 star HD 33959 (14 Aur) and the metal- line A7m star 
HD 88849 (HR 4021) found in The Elodie Archive at 
OHP. These stars are similar to DG Leo's components in 
spectral type and projected rotation velocity, but different in 
chemical composition. To make this more obvious, the spec- 
tra of HD 88849 and those of the Ab and B components were 
divided by the reference spectrum of HD 33959 and reported 
in Fig. [HI According to this comparison, DG Leo B is found 
to have almost the same chemical composition and funda- 
mental parameters as the reference star which is known to be 
chemically normal (|Hui-Bon-Hoall200fl ) . For their part, the 
A components present typical Am- type peculiarities (e.g.: 
scandium deficiency and enhanced iron-peak elements) but 
somewhat weaker than those observed in HD 88849. From 
a quantitative point of view, our abundance analysis clearly 
shows that in the A components' atmospheres, iron-peak el- 
ements are overabundant by up to 0.2 dex relatively to the 
Sun and by up to 1.0 dex for the heavier elements (Y, Zr, 
Nd). The overabundances are therefore of the same order 
of magnitude as is observed in other typical Am stars such 
as 32 Aqr (HD 209625, see Fig. [KJ. Deficient species, espe- 
cially calcium, are less extremely underabundant in the close 
binary components of DG Leo. 



8.5 The close binary components Aa and Ab 

Spin-orbit synchronization is generally expec ted to occu r 
100 times faster than orbital circularization jZahnlll977t) . 
Stars in an already circularized orbit are therefore expected 
to be synchronized and to have similar projected rotation ve- 
locities. As mentioned earlier, this is obviously not the case 



for the close binary of DG Leo. We note indeed (Sect. 6.1) 
a significant disagreement between the Vsini values of the 
two components. The predicted synchronization velocity 
V sy nch=36±6 km s _1 lies in between the measured values. 
Therefore, we cannot exclude that one of the components 
would be synchronized. This certainly merits further atten- 
tion. 

Components of multiple systems such as DG Leo usu- 
ally originate from the same protostellar environment. Their 
similarity or difference in chemical surface composition 
therefore relates directly to their individual/uncoupled evo- 
lution after formation (stellar evolution, rotation, diffusion), 
and the possible influence of a close companion. It is known 
from Sect. |S|that the Aa and Ab components have simi- 
lar fundamental parameters. As can be noted from Table |S] 
and from Fig. |K| they further have very similar abundance 
patterns. Both components are Am stars showing the same 
magnitude of overabundance for the iron-peak elements as 
well as for the heavier ones. The underabundances are mild 
compared to many classical Am stars. A differential com- 
parison of the disentangled spectra of Aa and Ab (Fig. |UJ 
gives evidence for a marginally higher Ca abundance in Aa. 
This fact remains unnoticed from Table |H| because the un- 
certainties given on the abundances include the uncertainty 
on the fundamental atmospheric parameters at an absolute 
level. When comparing Aa and Ab, only the relative errors 
in A T e g, Alog g and A£ tur b a re relevant. 

Classical diffusion models iAlecianll996t) generally show 
that calcium and scandi um are leaving the atmosphere after 
a while. Observations (lNorthlll993t iKunzli fc North! ll998T) 
and numerical simulations iAle*danlQl)96l) predict a short 
phase of calcium overabundance at the beginning of the 
main sequence followed by a calcium abundance decrease. 
In this sense, abundances are suspected to be very sensitive 
to mass loss, rotation and age which may explain the some- 
what different behaviour regarding diffusion observed in the 
components Aa and Ab of DG Leo. 



8.6 Spectral line variability 

The residuals of korel clearly indicate that the spectro- 
scopic variations due to pulsations are centred on the spec- 
trum of component B (Fig. 1101 . The pattern of these varia- 
tions is reproduced with a regularity of abou t 2 hours which 
agree s with the results from the photometry llLampens et alJ 
2005). A detailed study of the line variability awaits the 
analysis of the photometry as both sources should be com- 
bined in a periodicity analysis (our time base for spec- 
troscopy is very short) to enhance the probability of iden- 
tifying the pulsation modes. However, we note already the 
multi-periodic character of the line variability from the lack 
of exact repetition of the variability patterns in different 
nights. 

Despite this line variability, we further notice that the 
time-averaged spectrum of each component of DG Leo is 
fully consistent with the combined observed spectra and no 
numerical artifacts were detected. The use of korel to dis- 
entangle the component spectra from the high-resolution 
data may thus be considered as highly successful in this 
specific case. 
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resid 

Figure 10. The residuals provided by KOREL are plotted in a gray 
scaled figure for 4 consecutive nights (January 11 - 15, 2003). 
Each variation is centred on the spectrum of component B and 
repeats with a period of about 2 hours. 



codes. All three components appear to have the same funda- 
mental parameters (effective temperature and surface grav- 
ity) but DG Leo B was found to have a solar-like chemi- 
cal composition, while the two components belonging to the 
close spectroscopic binary show an abundance pattern typ- 
ical of Am stars. A different level of underabundance for 
scandium and more strongly for calcium is noted between 
the two close components (Aa and Ab). These differences 
probably reflect the sensitivity of the diffusion processes to 
issues such as rotation, mass loss and evolution. We further 
noticed that, while the spectroscopic orbit is circularized, 
the apparent rotation of both A components is not yet syn- 
chronized with the orbital motion. DG Leo is therefore a 
very interesting clue for understanding issues such as rota- 
tion synchronization and orbital circularization. 

Our high-time-resolution spectroscopic data confirm 
the existence of line-profile variations (LPVs) caused by pul- 
sation in component B with time scales ver y close to those 
recen tly detected by photometric studies llLampens et alJ 
2005). The analysis and the mode identification of these 
LPVs will be combined with the multicolour photometric 
observations, with the goal to provide a better understand- 
ing of the non-trivial link existing between chemical compo- 
sition, multiplicity and pulsation. 
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9 CONCLUSIONS 

The Fourie r Transform spe ctral disentangling technique de- 
veloped bv lHadraval lll995t) was successfully applied to the 
study of the pulsating triple system DG Leo. It allowed to 
confirm the pr eviously derived par ameters of the spectro- 
scopic orbit bv lFekel fc Bopd ( 1977) and to increase the ac- 
curacy. Combining astrometric measurements with our new 
radial-velocity measurements, we were further able to com- 
pute for the first time the relative orbit of the visual sys- 
tem AB. However, we stress the fact that speckle and RV 
measurements at periastron would be invaluable. We would 
like therefore to be informed about such data, if existing, 
as these measurements should greatly increase the precision 
of the present solution and it could further enable the de- 
termination of the stellar masses of all 3 components in an 
independent way. 

A detailed chemical analysis was performed on the dis- 
entangled spectra of the components with the Kurucz LTE 
model atmospheres combined with the SYNSPEC computer 
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